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Abstract 

We consider a model of hydrogen-like artificial molecule formed by two vertically coupled quantum dots in the 
shape of axially symmetrical thin layers with on-axis single donor impurity in each of them and with the magnetic 
field directed along the symmetry axis. We present numerical results for energies of some low-lying levels as 
functions of the magnetic field applied along the symmetry axis for different quantum dot heights, radii, and 
separations between them. The evolution of the Aharonov-Bohm oscillations of the energy levels with the increase 
of the separation between dots is analyzed. 
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Background 

An important feature in low-dimensional systems is 
the electron-electron interaction because it plays a crucial 
role in understanding the electrical transport properties of 
quantum dots (QDs) at low temperatures [1]. Such sys- 
tems may involve small or large numbers of electrons as 
well as being confined in one or more dimensions. The 
number of electrons in a QD can be varied over a consid- 
erable range. It is possible to control the size and the 
number of electrons and to observe their spatial distribu- 
tions in QDs. Energy spectrum of two-electron QD with a 
parabolic confinement, for which two-particle wave equa- 
tion can be separated completely, has been analyzed previ- 
ously by using different methods [2-5] . 

In the present work, we propose another exactly solv- 
able two-electron heterostructure in which two sepa- 
rated electrons are confined in vertically coupled QDs 
with a special lens-like morphology. Together with two 
on-axis donors, these two electrons generate an artificial 
hydrogen-like molecule whose properties can be con- 
trolled by varying the geometric parameters and the 
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strength of the magnetic field applied along the sym- 
metry axis. 

Methods 

The model which we analyze below consists of two 
identical, axially symmetrical and vertically coupled QDs 
with the on-axis donor located in each one of them (see 
Figure 1). The dimension of the heterostructure is defined 
by the QDs' radii R, height W, and the separation d be- 
tween them along the z-axis. We assume that the QDs 
have a shape of very thin layers whose profiles are given 
by the following dependence of the thickness of the layers 
w on the distance p from the axis: 



w(p) = W/Jl+(p/R) 



(1) 



Besides, for the sake of simplicity, we consider a model 
with infinite barrier confinement, which is defined in cy- 
lindrical coordinates as V(r) = 0 if 0 < z < w(p) , and 
V(r) = oo otherwise. 

Given that the thicknesses of the layers are much 
smaller than their lateral dimensions, one can take ad- 
vantage of the adiabatic approximation in order to ex- 
clude from consideration the rapid particle motions 
along the z-axis [6,7] and obtain the following 
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Figure 1 Scheme of the artificial hydrogen-like molecule. 



expression for effective Hamiltonian in polar coordi- 
nates: 
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The effective Bohr radius = h 2 e/m * e 2 as the unit 
of length, the effective Rydberg Ry* = e 2 /2eaQ = 
h 2 /2m * «g 2 as the energy unit, and y = ehB/2m * cRy* 
as the unit of the magnetic field strength have been used 
in Hamiltonian (Equation 2), with m* being the electron 
effective mass and c, the dielectric constant. The polar 
coordinates p k = (p k , dk) labeled by k = 1,2 correspond 
to the first and the second electrons, respectively. It is 
seen that for the selected particular profile given by 
Equation 1, the Hamiltonian (Equation 2) coincides with 
one which describes two particles in 2D quantum dot 
with parabolic confinement and renormalized inter- 
action. It is well known that such Hamiltonian may be 
separated by using the center of mass, R = (p 1 + P2V2, 
and the relative, p = p x — p 2 coordinates [8]: 
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The wave function is factorized into two parts, 
^(R : p) = <2>(R)0(p), describing the center of mass and 



the relative motions, respectively. Meanwhile, the total 
energy splits into two terms depending on two radial 
Ng, n p and two azimuthal Lr, l p quantum numbers: 



E(N R ,L R ;n p ,l p 



■ E s (Ng,Lg)+2E p (n p ,l p ) 
{2N R +L R )co + 2E p (n p ,l p ) 



(4) 



where the first term represents the well-known expression 
for the exact energy levels of a two-dimensional harmonic 
oscillator, labeled by the radial (Ng = 0, 1, 2, . . .) and 
azimuthal (L R = 0 ± 1 ± 2, . . .) quantum numbers 
for the center of mass motion and the relative motion 
energy 2E p (n p ,lp) must be found solving the following 
one-dimensional Schrodinger equation: 



-u"(p) + V(p)u(p) = E p (n p ,l p )u(p); V(p) 

= u 2 p 2 /<i+(l 2 -l/4)/p 2 

- 3/p - 4/vV + 4^2 (5) 



In our numerical, work the trigonometric sweep 
method [8] is used to solve this equation. 

Results and discussion 

Before the results are shown and discussed, it is useful 
to specify the labeling of quantum levels of the two- 
electron molecular complex. According to Equation 4, 
the energy levels E(N R .L R ; n p , l p j can be labeled by four 
symbols N R ,L R ; n p , l p . Even and odd l p correspond to the 
spin singlet and triplet states, respectively, consistent 
with the Pauli Exclusion Principle. 

We have performed numerical calculations of energy 
levels of complexes with radii R between 20 and 100 nm 
for different separations between layers. In all presented 
calculation results, the top thickness W is taken as 
0.4 nm. In order to highlight the role of the interplay be- 
tween the quantum size and correlation effects in the 
formation of the energy spectrum of our artificial system 
different from the natural hydrogen molecular complex, 
we have plotted in Figure 2 the potential curves E(d) = 
E(N Rl L R \n p ,lp) +2/d, similar to those of the hydrogen 
molecule in which the complex energies with the elec- 
trostatic repulsion between donors included as functions 
of the separation d between QDs are shown. Comparing 
them with the corresponding potential curves of the 
hydrogen molecule, one can to take into account that in 
analyzing the structure here, the electron motion in con- 
trast to hydrogen molecule is restricted inside two sepa- 
rated thin layers. The energy dependencies of different 
levels (labeled by four quantum numbers, (N Ry L R ; n p , l p j 
are shown in Figure 2 for QDs with two different radii, 
R = 40 nm and 7? = 100 nm. A clear difference in the be- 
havior of the potential curves is readily seen. If the 
curves are smooth without any crossovers for QDs of 
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Figure 2 Energies E(d) of the double-donor complex 
corresponding to some low-lying levels in vertically coupled 
QDs. As functions of the distance between them. 
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Figure 3 Energies E(N R , L R ; n p . I p ) of some low-lying levels of 
the double-donor complex in vertically coupled QDs. As functions 
of the magnetic field. 



small radius, the corresponding potential curves suffer a 
drastic change as the QD radius becomes large. In the 
last case, the energy levels become very sensitive to the 
variation of the separation between QDs, and the 
quantum size effect becomes essential, providing alter- 
ation of the energy gaps, multiple crossovers of levels 
with the same or different spins, and the level reorder- 
ing, as the distance between QDs increases from 5 to 
20 nm. 

We ascribe a dramatic alteration of the potential 
curves with the increase of the separation between QDs 
from 5 to 20 nm observed in Figure 2 to the interplay 
between the structural confinement and the electron- 
electron repulsion. As the QDs' radii are small(i? — > 0), 
the confinement is strong, and the kinetic energy (-1/ 
R 2 ) is larger than the electron-electron repulsion energy 
(~1/R), vice versa for QDs with large radii. Therefore, as 
the QDs' radii increase, the arrangement of the elec- 
tronic structure for different energy levels changes from 
typical for the gas-like system to crystal-like one, accom- 
panied by the crossovers of the curves and reordering of 
the levels. As the two-electron structure arrangement 
for large separation between electrons becomes almost 
rigid, the relative motion of electrons is frozen out, and 
the two-electron structure transforms into a rigid rotator 



with practically fixed separation between electrons. The 
electrons' motion in this case becomes similar to one in 
ID ring, and therefore, the energy dependencies on the 
external magnetic field applied along the symmetry axis 
should be similar to those which exhibit the Aharonov- 
Bohm effect. 
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Figure 4 Density of states for two different values of the 
magnetic field. Corresponding to low-lying levels of the double 
donor complex in vertically coupled QDs. 
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Figure 5 Density of states for three different distances between 
layers. Corresponding to low-lying levels of the double-donor 
complex in vertically coupled QDs. 



In order to verify this hypothesis, we present in Figure 3 
the calculated molecular complex energies E(Ng, Lr; n p , l p j 
of some lower levels as functions of the magnetic field 
strength for QDs with small R = 40 nm (upper curves) and 
large R = 100 nm radii (lower curves). 

It is seen that for QD of small radius, the energies are 
increased smoothly with very few intersections. Such de- 
pendence is typical for gas-like systems where the para- 
magnetic term contribution is depreciable in comparison 
with the diamagnetic one. On the contrary, the energy de- 
pendency curves for QD of large radius present multiple 
crossovers and level-ordering inversion as the magnetic 
field strength increases from 0 to 1. It is due to a competi- 
tion between diamagnetic (positive) and paramagnetic 
(negative) terms of the Hamiltonian whose contributions 
in total two-electron energy in QDs of large radii are of 
the same order while the electron arrangement is similar 
to a rigid rotator. In other words, the correlation in this 
case becomes as strong as the electrons are mainly located 
on the opposite sides within a narrow ring-like region. 

Finally, in Figures 4 and 5, we present results of the 
calculation of the density of electronic states for double- 
donor molecular complex confined in vertically coupled 
QDs. It is clear from the discussion above that the pres- 
ence of the magnetic field should provide a significant 
change of the density of the electronic states as the QDs' 
radii are sufficiently large. Indeed, it is seen from Figure 4 
that under relatively weak magnetic field (y = 0.5), as the 
molecular complex is confined in QDs of 100-nm with 
6-nm separation between them, the density of states 
becomes essentially more homogeneous since the widths 
of individual lines are broadened and the gaps between 
them are reduced. Such change of the density of states is 
observed due to a splitting and displacement of the indi- 
vidual lines accompanied by their crossovers and the 
reordering of the energy levels. 



In Figure 5, we present similar curves of the molecular 
complex density of states for three different separations 
between QDs. It is seen that the curves of the density of 
states are modified only slightly, essentially less than 
under variation of the magnetic field. Particularly, the 
lower energy peak positions are almost insensitive to any 
change of the distance between dots, while the upper en- 
ergy peaks are noticeably displaced toward higher energy 
regions. 

Conclusions 

In short, we propose a simple numerical procedure for 
calculating the energies and wave functions of a molecu- 
lar complex formed by two separated on-axis donors 
located at vertically coupled quantum dots with a par- 
ticular lens-type morphology which produces in-plane 
parabolic confinement. We show that in the adiabatic 
approximation, the Hamiltonian of this two-electron sys- 
tem included in the presence of the external magnetic 
field is separable. The curves of the energy dependencies 
on the external magnetic field and the separation be- 
tween quantum dots are presented. Analyzing the curves 
of the low-lying energies as functions of the magnetic 
field applied along the symmetry axis, we find that the 
two-electron configuration evolves from one similar to a 
rigid rotator to gas-like as the dot radii decrease. This 
quantum size effect is accompanied by a significant 
modification of the density of the energy states and the 
energy dependencies on the external magnetic field and 
geometric parameters of the structure. 
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